Phase coherence of an atomic Mott insulator 
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We investigate the phase coherence properties of ultracold Bose gases in optical lattices, with 
special emphasis on the Mott insulating phase. We show that phase coherence on short length 
scales persists even deep in the insulating phase, preserving a finite visibility of the interference 
pattern observed after free expansion. This behavior can be attributed to a coherent admixture of 
particle/hole pairs to the perfect Mott state for small but finite tunneling. In addition, small but 
reproducible "kinks" are seen in the visibility, in a broad range of atom numbers. We interpret them 
as signatures for density redistribution in the shell structure of the trapped Mott insulator. 
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A fundamental aspect of ultracold bosonic gases is 
their phase coherence. The existence of long-range phase 
coherence, inherent to the description of a Bose-Einstein 
condensate in terms of a coherent matter wave, was ex- 
perimentally demonstrated in interferometric |2, l| or 
spectroscopic 4] experiments. More recently, attention 
has been paid to fundamental mechanisms that may de- 
grade or even destroy long-range coherence, for exam- 
le thermal phase fluctuations in elongated condensates 
or the supernuid to Mott insulator (MI) tran- 
sition undergone in optical lattices |9|, |l(J, . 

For a Bose-Einstein condensate released from an opti- 
cal lattice, the density distribution after expansion shows 
a sharp interference pattern . In a perfect Mott Insu- 
lator, where atomic interactions pin the density to pre- 
cisely an integer number of atoms per site, phase co- 
herence is completely lost and no interference pattern 
is expected. The transition between these two limiting 
cases happens continuously as the lattice depth is in- 
creased. In the supernuid phase, a partial loss of long 
range coherence due to an increased quantum depletion 
has been observed for lattice depths below the MI tran- 
sition 0,0,^^. Conversely, in the insulating phase, 
numerical simulations [lMllull3 predict a residual inter- 
ference, although long-range coherence and superfluidity 
have vanished. 

In this Letter, we revisit this question of phase coher- 
ence focusing on the insulating phase. We observe that 
the interference pattern persists in the MI phase, and 
that its visibility decays rather slowly with increasing lat- 
tice depth. We explain this behavior as a manifestation of 
short-range coherence in the insulating phase, fundamen- 
tally due to a coherent admixture of particle/hole pairs 
to the ground state for large but finite lattice depths. In 
addition, we also observe reproducible "kinks" in the vis- 
ibility at well-defined lattice depths. We interpret them 
as signature of density redistribution in the shell struc- 
ture of a MI in an inhomogeneous potential, when regions 
with larger-than-unity filling form. Finally, the issue of 
adiabatic loading in the lattice is briefly discussed. 
In our experiment, a 87 Rb Bose-Einstein condensate 



is loaded into an optical lattice created by three or- 
thogonal pairs of counter-propagating laser beams (see 
flo| for more details). The superposition of the lattice 
beams, derived from a common source at a wavelength 
Xl = 850 nm, results in a simple cubic periodic potential 
with a lattice spacing d = = 425 nm. The lat- 

tice depth Vb is controlled by the laser intensities, and is 
measured here in units of the single-photon recoil energy, 
Er — h 2 /2m\ 2 L « h x 3.2 kHz, where m is the atomic 
mass. The optical lattice is ramped up in 160 ms, us- 
ing a smooth waveform that minimizes sudden changes 
at both ends of the ramp. After switching off the op- 
tical and magnetic potentials simultaneously and allow- 
ing for typically t = 10 — 22 ms of free expansion, stan- 
dard absorption imaging of the atom cloud yields a two- 
dimensional map of the density distribution (integrated 
along the probe line of sight). 

Four such images are shown in Fig. ^i-d, for various 
lattice depths. The density distribution of these expand- 
ing clouds can be expressed as [l£l LU| 
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In Eq. {I]), the interference pattern is described by 
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where the operator a \ creates an atom at site z, and where 
w is the Fourier transform of the Wannier function iu(r^). 
The Fourier relation (0) shows that long-range phase co- 
herence, i.e. a correlation function {d\dj) slowly varying 




FIG. 1: Absorption images of an ultracold Bose gas released 
from an optical lattice, for various lattice depths: a 8 i?R, b 
14 E Rl c 18 E R , and d 30 E R . 
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across the lattice, is necessary to observe a sharp diffrac- 
tion pattern as in Fig. QJi. However, above the MI tran- 
sition (Fig. CJ}, c an d d), the interference peaks evolve 
into a much broader, cross-like structure which weak- 
ens with increasing lattice depth. This slow modulation 
corresponds to short-range coherence, i.e. a correlation 
function (a\aj) whose range extends over a few sites only. 

To extract quantitative information from time-of- flight 
pictures as shown in Fig. [IJ Eq. (QJ suggests using the 
usual definition of the visibility of interference fringes, 
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In this work, we measure the maximum density n ma vat 
the first lateral peaks of the interference pattern |l9(, 
(i.e. at the center of the second Brillouin zone), whereas 
the minimum density n m i n is measured along a diagonal 
with the same distance from the central peak (see inset 
in Fig.EK)- I n this way, the Wannier envelope is the same 
for each term and cancels out in the division, yielding the 
contrast of S alone (hence the second equality in Eq. Q). 
Four pairs exist for a given absorption image, and their 
values are averaged to yield the visibility. In previous 
studies of the MI transition [HI Q, the sharpness of 
the interference pattern was characterized by the half- 
width of the central peak. Such a measure is possibly 
sensitive to systematic effects, such as optical saturation 
and mean field broadening. We expect our measure of 
contrast to be much less sensitive to these effects, since 
it is calculated in regions of the image where the density 
is lower. 

We present here measurements of the visibility as a 
function of lattice depth (typically in a range 6 — 30 £?r) 
at a given total atom number. Each value was obtained 
as the visibility averaged over approximately 10 indepen- 
dent images. Different atom numbers (hence different 
filling factors) were investigated, ranging from 6 x 10 4 
to 6 x 10 5 . Two illustrative sets of data are shown in 
Fig. O corresponding to approximately 5.9 x 10 5 atoms 
(black circles) and 3.6 x 10 5 atoms (grey circles). For 
lattice depths larger than 12.5 £^r, the system is in the 
insulating phase 10]. Yet, the visibility remains finite 
well above this point. For example, at a lattice depth 
of 15 £?r, the contrast is still around 30%, reducing to 
a few percent level only for a rather high lattice depth 
of 30 Er. We will now show that such a slow loss in 
visibility is expected in the ground state of the system. 

As shown in 0, the physics of ultracold atoms in an 
optical lattice can be described by the Bose-Hubbard 
hamiltonian, given by the sum of a tunneling term, 
Ht = —t^2^j^a\aj, plus an interaction term, Hint — 

(hi — 1). Here hi = a\ai is the on-site number 
operator, t is the tunneling matrix element, the nota- 
tion (i, j) restricts the sum to nearest neighbors only, 
and U is the on-site interaction energy [ll[. In the ex- 
periments, an additional, slowly varying potential V ex t(r) 




5 10 15 20 25 30 




0.15 



CD 

> 

CD 

Q 



5 10 15 20 25 30 
Lattice depth (Er) 

FIG. 2: (a) Visibility of the interference pattern produced 
by an ultracold cloud released from an optical lattice. The 
two sets of data shown correspond to 3.6 x 10 5 atoms (grey 
circles) and 5.9 x 10 5 atoms (black circles). The latter curve 
has been offset vertically for clarity. Arrows mark positions 
where "kinks" are visible, (b) Numerical derivative of the 
above curves. 



is also present, and favors the formation of a "wedding 
cake" structure of alternating MI and superfluid shells 
EJ [3 E3 •> which reflects the characteristic lobes delim- 
iting the MI phases in the phase diagram of the Bose- 
Hubbard model [ll|. 

To better understand the origin of a finite visibility, 
we consider a homogeneous system with filling factor no- 
In the limit of infinitely strong repulsion, U/t — > oo, the 
ground state is what we call a "perfect" Mott insulator, 
i.e. a uniform array of Fock states, |\I/)mi = EL l n o)i- 
This corresponds to a uniform S = no and zero visibility. 
To a good approximation, the actual ground state for 
a finite ratio U/t can be calculated by considering the 
tunneling term as a perturbation to the interaction term. 
To first order mt/U, this yields 
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The ground state thus acquires a small admixture of 
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FIG. 3: Visibility of the interference pattern versus U/zt, the 
characteristic ratio of interaction to kinetic energy. The data 
are identical to those shown in Fig. 0(5.9 x 10 5 , black circles, 
and 3.6 x 10 5 atoms, grey circles). The former curve has been 
offset vertically for clarity. The lines are fits to the data in 
the range 14 — 25 Er, assuming a power law behavior (see 
text). 



"particle- hole" pairs (i.e. an additional particle at one 
lattice site and a missing one in a neighboring site), which 
restores short-range coherence and a corresponding weak 
modulation in the momentum distribution, <S(k) ex no — 
2n (n +l)t(k)/U, where t(k) = -2t J2v= x , y ,z cos(k u d) is 
the tight-binding dispersion relation. The corresponding 
2D visibility (integrated along one direction) is 



(5) 



In Eq. (0, z = 6 is the number of nearest neighbors in a 
3D cubic lattice. 

To compare with the experiment, we show in Fig. 
the visibility against U/zt in a log- log plot. For lattice 
depths Vb > 14 Er (corresponding to U/zt > 8), the data 
matches the inverse law expected from Eq. (0. This has 
been verified by fitting the data in this range to a general 
power law A(U/zt) a (solid lines in Fig.EJ). We obtain an 
average exponent a = —0.98(7) in agreement with the 
prediction (see Fig. Efc). In Fig. ^3, the fitted prefactor 
is plotted as a function of atom number. Inspired by Eq. 
J5J), we compare it to 4(n + l)/3, where n is the average 
filling factor calculated at a lattice depth of 30 E^ using 
a mean- field approximation 0, [2^ . We find that this 
extrapolation of Eq. Q to our trapped system indeed 
yields the correct order of magnitude (see Fig. QJ)). We 
thus consider the agreement between our experimental 
results and the simple relations derived above as a con- 
clusive evidence for the presence of particle-hole pairs, 
characteristic of the ground state of the Bose-Hubbard 
hamiltonian. 

In addition to the smooth decay discussed above, the 
visibility shows small "kinks" at specific lattice depths 
(indicated by arrows in Fig. EH)- They are systemati- 
cally observed in our data, and their positions are re- 
producible. In the derivative plot (Fig. Et>), they ap- 
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FIG. 4: Exponent a (a) and prefactor A (b) extracted from 
a power law fit A(U / zt) a to the visibility data in Fig. 
plotted versus total atom number. The solid line indicates 
the expected exponent a = — 1. In (b), we also indicate the 
prefactor expected for uniform MI with filling factor no = 1 
(dashed line) and no = 2 (dotted line), as well as an extrap- 
olation for the average filling calculated at a lattice depth of 
30 E R (solid line). 



pear as narrow maxima on a smoother background. We 
obtained the kink positions by taking the middle point 
between two adjacent gaussian peaks with negative am- 
plitudes fitted to the data. The most prominent kink 
occurs on average for a lattice depth of 14.1(8) Er, with 
a statistical error indicated between parentheses. For the 
largest atom numbers (4.2 x 10 5 and 6 x 10 5 ), a similar 
but much weaker kink is also visible around 16.6(9) 
(see upper curves in Fig. EJ). These values are close to 
14.7 E R and 15.9 E R , the lattice depths where MI re- 
gions with filling factor no = 2 or 3 are respectively ex- 
pected to form for our parameters 22]. We thus propose 
that the observed kinks are linked to a redistribution in 
the density as the superfluid shells transform into MI 
regions with several atoms per site. We were recently 
informed that similar features were reproduced numeri- 
cally for one-dimensional trapped systems with a small 
number of particles |23j |. 

We have considered the dependence of the visibility 
on the time over which the optical lattice was ramped 
from zero to its final value, for a specific lattice depth 
of Vb = 10 £?r. The visibility was considerably de- 
graded for the shortest ramp time of 20 ms, but reached 
a ramp-independent value for ramp times larger than 
T&d ~ 100 ms (to be compared to the 160 ms time used 
in visibility experiments). We note that T a( j for this lat- 
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tice depth of Vq = 10 Er is significantly longer than the 
microscopic time scales of the system, such as the tun- 
neling time or the trapping periods. We note also that 
at the largest lattice depth we use here (Vq = 30 £^r), 
the observed visibility is systematically above the power 
law fit in Fig. indicating a breakdown of adiabaticity. 
By comparing the data to the fitted curve, we expect 
this to occur for Vq « 29 Er (U/zt& 200), which agrees 
with the calculated depth of 32 for which the ramping 
time 160 ms becomes smaller than the calculated tunnel- 
ing time h/zt. 

Although a complete study is beyond the scope of this 
Letter, these observations suggest that different dynam- 
ical processes are involved in the loading, depending on 
whether the gas is in the superfluid or in the MI phase. 
In the superfluid phase, the ramp time has to be slow 
enough not to excite long-lived collective excitations. In 
the MI phase, these excitations acquire an energy gap, 
which makes single particle tunneling the dominant dy- 
namical process. In this case, the final tunneling time in- 
creases with final lattice depth, and eventually becomes 
so long that the system basically freezes out at some lat- 
tice depth, estimated here to be 29 

In conclusion, we have studied the visibility of the in- 
terference pattern produced by an ultracold Bose gas re- 
leased from a deep optical lattice. A non- vanishing visi- 
bility in the MI phase is observed and explained by the 
coherent admixture of particle-hole pairs to the insulat- 
ing ground state, which preserves local phase coherence. 
This intrinsic limitation to the "quality" of a MI has 
important implications for various quantum information 
processing schemes, where the MI plays a central role 
13 EH E3- In addition, we observe small but repro- 
ducible kinks in the visibility curve. We interpret them as 
the signature of density redistribution in the shell struc- 
ture of the cloud as MI with several atoms per site are ex- 
pected to form. Finally, a recent paper 27] suggests that 
in a planar array of one-dimensional Bose gases, the vis- 
ibility might be further reduced when correlations build 
up in each tube, i.e. upon entering the Tonks- Girardeau 
regime. Experimental study of these effects seems within 
reach with the methods presented in this paper. 
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